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jH . Abstract Aim of this paper is the qualitative analysis of a boundary value problem for a third 

^ . order non linear parabolic equation which describes several dissipative models. When the source 

term is linear, the problem is explictly solved by means of a Fourier series with properties of rapid 
convergence. In the non linear case, appropriate estimates of this series allow to deduce the asymp- 
totic behaviour of the solution. 
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^- 1. INTRODUCTION 

^ ■ The evolution of several dissipative phenomena can be described by the non linear equation 

Ceu = {d^^iedt + c2) - dt{dt + a))u ^ F{x, t, u), (1.1) 

where s,a,c are positive constants and F = F{x,t,u) is a prefixed function. For instance, the 

K^ ■ equation (1.1) governs sound propagation in viscous gases as reported in the classical literature 

^ I [1] , motions of viscoelastic fluids or solids [2] , or motions of heat condution at low temperature 

-^ ■ (see, e.g., [3] and refecence therein). Moreover, in superconductivity, (1.1) models the flux 

pg I dynamics in the Josephson junction. In this case, let u — u{x,t) be the phase difference of the 

wave functions related to the two superconductors, and let 7 be the normalized current bias; when 

^—v I F = sinu — 7, equation (1.1) is the perturbed Sine - Gordon equation (PSGE). Terms su^xt a-nd 

fSj ■ aut characterize the dissipative normal electron current flow along and across the junction and 

represent the perturbations with respect to the classic Sine Gordon equation [4]. Moreover, the 

value's range for a and e depends on the real junction. In fact, there are cases with < a, e < 1 

[5] and, when the resistence of the junction is so low to short completely the capacitance, the case 

a large with respect to 1 arises [6]. As it is well known, up to today, Josephson junctions have 

jrt i found a lot of applications in many fields and lately new high- Tc materials imply an increasing 

" " evolution of the superconducting technologies [7] . 

If i is the normalized length of the junction, boundary conditions for a: = and x = £ specify 
the applied external magnetic field [8, 9]. Then, if 

n^{{x,t) -.0 <x <i, 0<t<T}, 

the strip problem for (1.1) can be stated as follows: 

dxx{£Ut + cPu) ~ dt{ut + au) = F{x,t,u), [x,t) e SI, 

u{x,0) = goix), ut(x,0)=5i(x), xe[OJ], (1.2) 

u(0,t) = 0, u{i,t)^0, 0<t<T. 



The linear case(i^ = f{x, t)) has aheady solved in [10, 11] by means of convolutions of Bessel 
functions. However, it seems difficult to deduce an exhaustive asymptotic analysis by means of 
this solution. 

In this paper, both linear and non linear problems are analyzed and the Green function of 
the linear case is determined by Fourier series: 

r, OO 

G{x,t,£,) = -^Hn{t) sin7„^ sin7„x, (1.3) 



^.1 



where 



i7„(t) = — e-''"*sinhKi), -fn^"^, (1.4) 

UJn I 

bn = C7„, hn = -{a + e-fl), uj^ = ^hl - b^. 

This series is characterized by properties of rapid convergence and allows to establish an 
exponential decrease of the solution in the linear case. Moreover, the non linear problem is 
reduced to an integral equation with kernel G and, by means of suitable properties of G, the 
asymptotic analysis of the solutions is achieved. 

2. PROPERTIES OF THE GREEN FUNCTION 

Firstly, let us prove that function G defined by series (1.3) decreases at least with n~^ and 
that it is exponentially vanishing as i — > oo. In fact, if we put: 

P= e + aHl.f ' '^ = 2 ' /5 = mm(p,g), (2.1) 

the following lemma holds: 

Lemma 2.1- The function G{x,^,t) defined in (1.3) and all its time derivatives are continuous 
functions in fl, V a^EjC^ G 3?"'". Moreover, everywhere in fl, it results: 

|G(x,C,t)| < Me-^ \—-\<N,e-^\ j eU (2.2) 

where M,Nj are constants depending on a,e,c^. 

Proof. When c^ — ae, the operator L^ can be reduced to the wave operator. If c^ < ae, for 
all n > 1, one has ft-„ > 5„. So, if k is an arbitrary constant such that c^ /ae < fc < 1, Cauchy 
inequality assures that: 

^ <l{e^l + a) = h^. (2.3) 



Then, if X„ = (5„//i„)^, it results: 



{l-k)i <^ = {1-X,,)i <l-Xn/2. (2.4) 

hn 

So, considering that: 

^> P, hn>n^{q-a/2), Vn > 1, (2.5) 

g-t(/i„-w„) ^ ^~pt g^^^ 1^^ <- ( ^_\^ -^ hold. As consequence, (2.2)i follows. 

As for (2.2)2, we observe that hn^ijJn < -^, Vn > 1, and by means of standard computations 
(2.2)2 is deduced. 



It may be similarly proved that estimates (2.2) hold also when c^ > ae. I 

As for the x-differentiation of Fourier's series like (1.3), attention must be given to convergence 

problems. Therefore, we consider x-derivatives of the operator {edt + c?)G instead of G and G*. 

Lemma 2.2.- For all a^e^c^ G 3fJ+, the function G{x,(^,t) defined in (1.3), is such that: 

leGt+c'Gl <Aoe-^\ \d^^{eGt + c^ G)\ < Ai er^\ (2.6) 

where Ai (i = 0, 1) are constants depending on a,e,c^. 
Proof. As for the hyperbolic terms in iJ„, it results: 

£ff„ + c^H^ = -Le-'*"*{(c' + e<^n)e""* - [c" - e(hn + w„)]e-""*}, (2.7) 

with ipn = hn{-l + VI - Xn) and X„ = (^) < 1. 

Further, by means of Taylor's formula it's possible to prove that: 

\c^+£ipn\ < ^ (5aV4 + ai/8), (a,ai = const). (2.8) 

Estimates of Lemma 2.1 together with (2.8) show that the terms of the series related to the 
operator (2.6)i decrease at least as n^^. So, it can be differentiated term by terms with respect 
to X and estimate(2.6)2 holds. I 

Wc mean as solution of the equation CgV = a continuous function v{x, t) which has contin- 
uous the derivatives Vt,vtt, dxx{£Vt + c^v) and these derivatives verify the equation. Therefore, 
the following theorem holds: 

Theorem 2.1- The function G — G{x, t) defined in (1.3)is a solution of the equation 

CeG - dxAeGt + c^G) - dt{Gt + aG) = 0. (2.9) 

I 
To obtain the explicit solution of (1.2), properties for the convolution of function G with data 
have to be analized. 
For this, let g{x) be a continuous function on {0,C.) and consider: 

Ug{x,t) ^ g{^) G{x,^,t) d^, u*g{x,t) = {dt + a-edxx)ug{x,t). (2.10) 

Then one has: 

Lemma 2.3- If g(x) is a G^[{),1], then the function Ug is a solution of the equation L^ ~ 
such that: 

]iTaUg{x,t) ~ Q, lmidtUg{x,t) — g{x), (2-11) 

uniformly for all a; € [0,i]. 

Proof. Lemma (2.2)-(2.3) and continuity of g assure that function (2.10)i and partial deriva- 
tives required by the solution converge absolutely for all {x,t) G f2. Hence, since Theorem 2.1, 
CeUg = 0. Besides, hypotheses on the function g and (2.2), imply (2.11)i. Further, one has: 

2 a ,^ ■ cos7„^ . 

Gt = ^ >^ H„[t) sm7„x, (2.12) 

TT ot, _ n 



and hence: 



2 °° ff (f) 
dtUg = V" —^ — [5(0cos7„^]o sin7„x (2.13) 



TT 

n=l 



nl OO 



So, if 77(0;) is the Heaviside function, from (2.13), one obtains: 



lirndtug = ^[gie) - 5(0)] + g(0) - j^ 9 {Olvi^ " "^^ + i " ^^"^^ = ^^'^^^ ^^"^^^ 

I 
Lemma 2.4- //.g(a;) G C^p,^] wif/i g**'(0) = g^'\i) = {i ^ 1,2,3), t/ie^ f/ie function u* 
defined in (2.10)2 "is a solution of the equation C^ = such that: 

\im u*g{x,t) — g{x), \iindtu*{x,t) — 0, (2-15) 

uniformly for all x G [0, (]. 

Proof. Hypoteses on g{x) assure that: 



d^^ug{x,t)^ / g (0 G(x,C,t) di^u„{x,t). (2.16) 

JO 

So, by Lemma 2.3 and expression (2.16), equation C^u* = is verified. Moreover, since dtu* = 
c^Ugn and (2.11)i, (2.15)2 holds. Finally, owing to (2.11)-(2.16), one has (2.15)i, too. I 

3. EXPLICIT SOLUTION OF THE LINEAR PROBLEM 

Consider the linear problem. When F = , by Lemma 2.3 and 2.4, one has: 

Theorem 3.1- If the initial data gi{x), and go{x) verify hypotheses of Lemma 2.3 and 2.4, 
respectively, then the homogeneous problem (1.2) admittes the following solution: 

uo{x,t) ^Ug^ + [dt +a-ed^^)ug„. (3.1) 

I 
As for the linear non-homogeneous problem, (F = f{x,t)), let: 

Uf{x,t) =- f dr f /(e,T) G{x,U-r) d^. (3.2) 

Jo Jo 

By means of standard computations, one can verify that the function (3.2) satisfies (1.2) with 
go — gi — (^- As consequence, one has: 

Theorem 3.2- // / and f^ are continuous functions in 51, then the function u = Uq + Uf 
represents a solution of the linear non-homogeneous strip problem. 
Proof. Considering that: 

lim / f{i,T)Gt{xX.t-T)d^^f{x,t), (3.3) 

Theorem 2.1 assure that C^Uf = f(x,t). 
Moreover, since: 

|u/|<5i(l-e-'5*); \dtUf\<B2{l-e-^'), (3^,82 -const.), (3.4) 

initial homogeneous conditions are satisfy. I 

Uniqueness is an obviously consequence of the energy- method (see, for example, [12]). So, 
the following theorem holds: : 



Theorem 3.3- When the source term f{x,t) satisfies Theorem 3.2, and the initial data 
{go,gi) satisfy Theorem 3.1, then the function 

u{x, t) = Ug^ + {dt+ a- sdxx)ug„ + Uf (3.5) 

is the unique solution of the linear non-homogeneous strip problem with F = f{x,t). I 

Consider, now, the non-linear piohlcia (1.2), and observe that (3.5) imphes: 

u{x,t)^ f gi{OG{x,tt)d^ + {dt + a-edxx) f goiOGix,U)d^ (3-6) 

Jo Jo 

dr / G{x,tt-T)F{^,T,u{^,T))d^. 
Jo 

which represents an integral equation for u(x,t). 

4. ASYMPTOTIC PROPERTIES 

Obviously, asymptotic properties depend on the shape of the source term F. Therefore, when 
F = 0, Lemma 3.1 assures that the solution exponentially vanishes when t tends to infinity. 
When F — f{x,t) one has: 

Theorem 4.1- When the source term f{x,t) satisfies condition: 



1 

[k + ty 



I /(x, t)\<h ,,^ , {h, k,a^ const > 0), (4.1) 



then solution (3.5) is vanishing as t —> oo, at least as t ". 
Moreover, when 

\f{x,t)\ < Ce-^\ {C,S = const >0), (4.2) 

one has: 

\u{x,t)\ < ke^^"*, S* ^min{(3,6}, fc = const. (4.3) 

Proof. Since Ug^ ~ 0{e^^^) (i=0,l), the proof of the theorem follows directly from the 
assumptions on /(x,t). I 

As for the non linear problem, previous estimates can be applied. For instance, if we refer to 
PSGE, it results: 

|i^(a;,i,w)| = I sinu + 7I < 71, 71 = const. (4.4) 

Then, the solution of the related non linear strip problem is bounded for all t. 

Results concerning the exponential decay can be obtained when stronger assumptions on the 
decreasing properties of the source F(x,t,u(x,t)) are imposed. For instance, according to [13], 
when 

|F(x,t,u(x,i))| < const, e"''* (m > 0), (4.5) 

then, the solution u(x,t) of the non feear problem (1.2) vanishes exponentially, due to properties 
(2.2) of G. 
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